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ABSTRACT 

A software program has been developed and tested to process 
microwave radiometry data to be acquired by the microwave sen- 
sor (SSM/T) on the Defense Meteorological Satellite Program 
(DMSP) spacecraft. The SSM/T 7-channel microwave radiom- 
eter and systems data will be data-linked to Air Force Global 
Weather Central (AFGWC) where they will be merged with ephem- 
eris data prior to product processing for use in the AFGWC 
Upper Air Data Base (UADB). The overall system utilizes an 
integrated design to provide atmospheric temperature soundings 
for global applications. The fully automated processing at AFGWC 
is accomplished by four related computer processor programs 
to produce compatible UADB soundings, evaluate system per- 
formance, and update the a priori developed inversion matrices. 
Tests with simulated data produced results significantly better 
than climatology. 


1. INTRODUCTION 

With the SSM/T temperature sounder, the United States Air Force Space and Missile 
Systems Organization (SAMSO) and AFGWC are continuing efforts to improve the accuracy of 
global atmospheric temperature soundings. The objective is to overcome the limitations of 
infrared radiometers by providing a better capability to recover such soundings xinder cloudy 
conditions. This is expected to contribute significantly to AFGWC weather forecasting 
capability. 

The SSM/T software is designed to satisfy the specific requirements of AFGWC. They 
include requirements for the operating environment, inputs, data location, operator interac- 
tion, calibration data correction, inversion algorithm, and output products. In the opera- 
tional software, the brightness temperatures are processed to yield atmospheric tempera- 
tures for mandatory pressure levels of 1000, 850, 700, 500, 400, 300, 250, 200, 150, 100. 
70, 50, 30, 20 and 10 mb; 14 thicknesses between the mandatory levels; and the temperature 
and pressure of the tropopause. These comprise the parameter vector defined in Section 3. 
The heights of the mandatory levels are computed by stacking the thicknesses on to the fore- 
cast 1000-mb height which is retrieved from the AFGWC system. These parameter data are 
stored in catalogued disk files prior to data transfer and archiving for diagnostic analysis 
and updating purposes. 

To obtain these parameters, the SSM/T sounder system, developed at Aerojet Electro- 
Systems Company, operates in the oxygen absorption band at frequencies of 50. 5, 53. 2, 

54. 35, 54. 9, 58. 4, 58. 825 and 59. 4 GHz. The sensor is a crosstrack scanning radiometer 
which acquires data at 32-second intervals at 7 angular positions separated by 12 degrees. A 
multiple regression technique, applied after the radiation temperatures are adjusted for 
antenna pattern and terrain height effects, is used to invert the microwave data. The inver- 
sion matrices (one for each season, geographic zone and scan position) have been calculated 
from carefully selected a priori ensembles of atmospheric soundings and preflight-calibrated 
radiometer characteristics. Provisions are made to update and correct the inversion matrices 
cifter launch using coincident microwave and composite radiosonde/rocketsonde observations. 

The SSM/T sensor system is expected to be operational in 1978. In the meantime, the 
software system has been developed, tested and implemented into the AFGWC system at 
Offutt Air Force Base, Nebraska. The tests, involving the use of simulated data, indicated 
that soundings will be obtained with an accuracy considerably better than estimates based on 
climatology. 
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2. SYSTEM DESCRIPTION 


The SSM/T sensor system is designed to acquire microwave data simultaneously in the 
seven channels shown in Table I. Design specifications for these channels are also given in 
this table. Channel 1 is a window channel responding strongly to the earth's surface charac- 
teristics, dense clouds, and rain. It is centered at a frequency of 50. 5 GHz and is used to 
correct the other channels for these background effects. Since the polarization dependent sur- 
face effects are removed from the brightness temperature measurements, only the atmos- 
pheric effects, which are the same for all polarizations, remain. Thus, a priori polarization 
corrections are avoided. To achieve this result. Channels 1-4 are required to receive the 
same, but arbitrary, polarization. Since Channels 5-7 are not affected by the surface, their 
polarization may be arbitrary. The choice of an orthogonal mode was made solely for the 
sensor frequency separation design requirements. Frequency stabilization, stepped auto- 
matic gain control (SAGC), and instrumentation to provide critical component temperature 
data are further refinements of the sensor design that have a bearing on the software. 

Microwave data are acquired by means of a mechanically scanning, shrouded-reflector 
antenna system. The calibration and data acquisition scanning geometry are illustrated in 
Figure 1 and the scan parameters in Table II. A planar scan is generated by mechanically 
rotating the shrouded reflector perpendicular to the flight path. As shown in Figure 1, micro- 
wave data are acquired for all frequencies at seven scan positions at 12 -degree angular inter- 
vals. On each scan the radiometer is calibrated; first by a 300 K warm reference load and 
then by sensing the cold space cosmic background (2. 7 K). The scan period is 32 seconds 
when synchronized and 29. 4 seconds when in automatic mode. Synchronization, control, 
clock reference, system time and source power are furnished by the DMSP spacecraft. 

For each frequency, the antenna temperature, T^ , is determined in the software by the 
fundamental relationship 


■^r - 


H 


( 1 ) 


where 


Tj^ = warm calibration reference load input temperature 
T^ = cold space calibration input temperature 

= output voltage during warm load calibration 
= output voltage during cold space calibration 
V = output voltage corresponding to T^ 


3. RETRIEVAL TECHNIQUE 


The SSM/T retrieval method is based on a linear multiple regression technique. This 
method, developed specifically for the DMSP SSM/T application, is described briefly below. 

It is assumed that the deviation of the parameter vector from the climatological mean may be 
expressed as a linear combination of the deviations of the measured data from its mean. Thus, 


where 



D(d-<d>) 

atmospheric parameter vector, d = data vector 
predicted value of p , = expected value of p 

expected value of the data vector, and D is a matrix. 


( 2 ) 


Equation (2) requires the matrix D to have as many rows as p has components and as many 
columns as d has components. The matrix D is determined by the condition that the ex- 
pected value of the square of the difference between the true value of the parameter vector p 
and p be a minimum. This leads to the equation 


D = C|p-<p>, d -<d>)c"^(d -<d), d-<d>) 

where C(. , . ) denotes the covariance matrix of its two arguments. 


( 3 ) 
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where 


By collecting terms, Equation (2) may be rewritten in the form 
p = Dd + A 
A = <p> - D<d> 


(4) 

(5) 


A and a linear combination 


so that the predicted parameter vector is expressed as a constant 
of the data. 

To determine the covariance matrices required in Equation (3), either real time data 
gathered during sensor operation or a priori data may be used. Since a large number of upper 
air soundings, representing a wide range of meteorological conditions, have been archived 
over the years and the microwave response to a given atmosphere may be calculated, the 
latter method was considered to be more practical for implementation. 

Since microwave radiometers possess inherent instrumental noise characteristics, the 
noise must be accounted for in the computations. This is accomplished by writing the data 
vector d as 


d = d + 8 
o 

where d^ is the exact value and 6 a random error representing system noise with an 
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expected value of zero. By replacing 
covariance matrices are found to be 


d on the right side of Equation (3) by Equation (6), the 


C (p-<P>, d-<d>) = C (p - <p) , d^ - 

c(d-<d>,d-<d>) = c(d^-^d^), d^-^d^))+C(6, 
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Equation (8) incorporates the expected instrumental noise characteristics in the computations 
through the noise covariance matrix C(6, 6) . 

As indicated above, the data vector d required for use in the multiple regression method 
is constructed from brightness temperatures. Frequencies at which the brightness tempera- 
tures are measured are chosen to respond to different atmospheric layers ranging from sea 
level to the 10 mb level. The resulting vertical coverage is determined by weighting functions 
which indicate the amount of microwave radiation arising from specific altitude regimes. 

The weighting functions are defined by the following procedure. The radiative transfer 
equation shows that the brightness temperature may be expressed as 


T = L 
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g sky^ 
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e = 


g 


angle with respect to earth's surface 

air temperature at height Z 

absorption per unit length in the atmosphere 

emissivity of the ground 

ground thermal temperature 


Equations (10)-(12) are the transmission factor, the downward flowing component of the bright- 
ness temperature at the surface, and the upwelling component emitted between the earth and 
the microwave sensor. The contribution of the radiation emitted by the atmosphere to the 
temperature T is isolated from the ground component term by writing Equation (9) as 
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T = L eT + 
g 


T . dZ 
air 


where 'N(2.) is the weighting function given by 


L sec0 « <( 1 - e ) exp 


r r ^ 1 

-secS j* «dZ + exp sec^ J" »dZ 
o J L o . 


The function W specifies the contribution of the atmosphere in a thickness dZ at an 
altitude Z to the total signal energy received. Weighting functions corresponding to the 7 
SSM/T frequencies are shown in Figure 2. These are computed for an ARDC model atmos- 
phere and a calm sea background corresponding to an incident angle of 0° using a nominal 
spacecraft orbit altitude of 450 nmi. 

For the DMSP application a data vector was chosen to depend only on the atmospheric 
state. This data vector is derived from a rearrangement of Equation (9) as 


LT , +T = T - «T L 
sky atm g 




where 


d. = (L T , + T ^ ) 

i sky atm i/j 

= radiometer frequencies. 


Since the right-hand side of Eouation (16) is not directly measurable, it must be com- 
puted from the measurements. This is made possible by the proper employment of the 50. 5 
GHz channel to determine *Tg and the lower troposphere contributions from dense clouds 
and rain. Thus the data vector defined in Equation (16) is a 6-component vector where i/j 
are the 6 highest frequencies previously specified. 

To fully understand the ramifications of the SSM/T data vector concept. Equation (9) is 
rewritten as 


<‘V50 


. 5 = [^ - 


L T , - T ^ 

sky atm 




specialized to 50, 5 GHz. Since the same portion of the earth is viewed by all radiometer 
channels and with the same polarization. Equation (17) is substituted in Equation (15) to pro- 
vide a measured data vector defined as 


di Ty - (T^q g - s^Q g) a^^ 

I I 

®50. 5 ^ ^sky ^atm^50. 5 


(i = 1. 2. 




(i = 1. 2. . . 


These expected values are computed from the same a priori data used to determine the 
covariance matrices. 

When these equations are applied to the form of the multiple regression method described 
by Equations (4) and (5), the estimate of the parameter vector p is given by 

p = S Tg + a' (19) 

where S = ( -Da! D) is a matrix whose first column is -Da and whose remaining columns are 
the colunrms of D , 
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and 


A' = <P> - D<d>+ S 50 5 Da 

In summary this retrieval technique contains elements which depend mainly on the atmos- 
phere not on the background. It is equally valid over land, water or mixed surface conditions. 
In contrast to the ground, the atmosphere emits unpolarized radiation so that the received 
polarization is no longer a factor in the data processing. The data processing does provide, 
however, corrections for antenna pattern and terrain height effects. These are described below. 

Antenna pattern corrections are applied to the antenna temperatures directly after cali- 
bration to yield brightness temperatures. The method consists of scaling theantenna tempera- 
ture by a factor gg which is the fraction of the total gain of the anPenna intercept- 

ing the earth, to produce the brightness temperature. 

A total of 49 factors are used, one for each frequency and scan angle. 

Since the inversion matrices are computed only for sea level surface conditions, adjust- 
ments to the brightness temperatures determined above must be made for measurements 
taken over elevated terrain. This adjustment is necessary because the brightness tempera- 
tures measured by the lower microwave frequency channels may differ significantly from 
temperatures which would arise were the same surface at sea level. The effect caused by the 
missing portion of the atmosphere requires an adjustment to the brightness temperatures 
prior to inverting. The adjustment is derived from Equations (9) and (17) as 

T (u) = T ( 1 ^) + A - €T A, 

o 1 a g b 

where 

A(i') = <T, +LT,\-/t. +LT,\, 

a \ atm sky/o \ atm sky/ 1 

and 

In these equations, the subscript 1 represents the true atmosphere and surface at height 
h while o represents the true atmosphere above h and an intervening atmosphere which is 
determined statistically between sea level and h . The quantities 850 5 and the first term 
in the denominator of Equation (24) have been defined in Equations (18) and (17) respectively. 
The factors Aa(t') and Ai,(i/) are computed only for Channel 1 to 4 frequencies since the 
three highest frequency channels are unaffected by the terrain. These factors are determined 
by the characteristics of the atmosphere comprising the statistical data ensembles. Quad- 
ratic functions were found to fit the data satisfactorily. Specifically, for each of the low 
channels 

Aa = A^h + Aj^h^ (25) 

and ~ 

A. = B h + B, h'^ 
b a b 

where h is the terrain height and the constants A^ , A{j , B^ , Bj, are regression values 
determined from the data base ensembles. 

The SSM/T software contains provisions to correct the a priori D-matrices during the 
operational phase to optimize the retrieved data when necessary. The procedure consists of 
comparing SSM/T derived soundings with coincident radiosonde soundings. The implemented 
update procedure corrects for sensor and atmospheric modeling bias errors. A linear 
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(23) 
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+ Aj^(50. 5) 


(24) 
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regression is performed between the measured data vector and a data vector calculated from 
coincident composite radiosonde and rocketsonde observations. Regression coefficients 
and are computed for the ith frequency by 

d = ad + /3 (27) 

me 


where dj^j is the measured data vector and dj. the calculated data vector. These coeffi- 
cients are used to update the a priori correlation matrices (indicated by the subscript a) and 
the expected data vector <da_^ to yield updated quantities (indicated by the subscript u) as 
shown in Equations (28) -(31) 


P = 




C(p-<P>, = C (p - <^p> . d^ - <d^>) a 

c(V<'*u>'V<0) - '*a-<0)“ 

<du> = «<d^>+ P 


(28) 

( 29 ) 

(30) 

(31) 


These updated matrices are combined with an operationally determined error covariance 
matrix, C(S,8) , to yield new matrices defined by 


D = C(p -<p> , d^ - <d^) [c(d^ - <dj>, d^ - <d^>)+ C(6,6)]-' 


(32) 


4. DATA PROCESSING 


4. 1 GENERAL DESCRIPTION 

Ephemeris data are merged with the raw microwave sensor data and DMSP system data 
at AFGWe to form the required input files for processing. The SSM/T data processing con- 
sists of four main computer program components labeled PROCESSOR-A, PROCESSOR-B, 
MICRO-CHECK, and D-UPDATE. Only PROCESSOR-A and PROCESSOR-B are designed to 
operate in real time. PROCESSOR-A is fundamentally a preprocessor for PROCESSOR-B, 
the inversion processor. MICRO-CHECK is an on-line program, operating once each day as 
a preprocessor for the update program, D-UPDATE. Both MICRO-CHECK and D-UPDATE 
are designed to evaluate and correct the a priori D-matrices. Each of these processing pro- 
grams consists of a driver routine and numerous subroutines. They are described in the 
following section in reference to the top level flow diagram shown in Figure 3. 

4.2 PROCESSOR-A 

PROCESSOR-A is designed to process data contained in the input file PROD^AEROJET- 
MDATA while retrieving required information from various AFGWC data base components. 
The processing is initiated at the UNIVAC 1110 system console by operator keyboard input. 
Data are read in convenient groups from the input file to minimize input/output operations. 
Raw data, in digital counts, are processed to produce the SSMT*BRTEMP file, SSMT*- 
CALGAIN file, and various page print diagnostics. SSMT*BRTEMP file serves as the main 
input to PROCESSOR-B whereas SSMT*CALGAIN is generated for off-line operational analy- 
sis and system evaluation. The diagnostics are produced for error analysis purposes to aid 
in the recovery of possible program failures. As depicted in Figure 3, valid sensor data are 
stripped from the input file and identified for each spacecraft. These data are calibrated, 
corrected for antenna patterns, timed, and earth-located. Gross error checks are made to 
the corrected brightness temperatures against fixed upper and lower limits. Corresponding 
terrain heights and 1000-mb heights are retrieved from appropriate elements in the AFGWC 
data base. The corrected brightness temperatures and associated view data are stored in 
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the SSMT*BRTEMP file while validated raw data are stored in the SSMT*CALGAIN file. 
Page-print diagnostics may be generated for this processor by the operator. PROCESSOR -A 
executes automatically processing one or more orbits of data from the input file which may 
contain a maximum of 24 hours of data. After PROCESSOR -A executes, further real-time 
processing is continued in PROCESSOR-B 

4.3 PROCESSOR-B 

PROCESSOR-B runs sequentially to PROCESSOR -A in an automatic mode with the proc- 
essing mode and area limits keyed in by the operator. Processing involves reading data from 
the SSMT*BRTEMP file, selecting data for appropriate views (scan position), testing for the 
likelihood of dense clouds or rain (over oceans only), adjusting for terrain height and invert- 
ing to obtain atmospheric temperature soundings. The inversion process is described in 
Section 3. All a priori D-matrices are stored in the SSMT*D-MATR1CES file. There are 48 
D-matrices in all, including one each for three geographic zones (tropics, mid-latitudes and 
arctic /antarctic), four seasons, and four scan positions. The recovered soundings are stored 
in the SSMT*SOUN DINGS file and SSMT^VERIFICATION file. The corresponding brightness 
temperatures are also stored in the SSMT^VERIFICATION file. In addition PROCESSOR-B 
produces page-print diagnostics of files as well as the location of views having dense clouds 
or rain. After processing the last view of the SSMT*BRTEMP file, PROCESSOR-B provides 
a listing of the number of soundings processed prior to termination. 

4. 4 MICRO-CHECK 

The MICRO-CHECK program is designed to run once each day. Inputs to MICRO-CHECK 
are derived from the SSMT*VERIFICATION and AFGWC UADB. MICRO-CHECK searches 
the UADB for radiosonde observations coincident with the inverted SSM/T soundings. It com- 
putes the residual errors between the SSM/T soundings and the coincident radiosonde obser- 
vation. The residual error data are used to evaluate the overall performance of the retrieved 
profiles. Coincident rocketsonde soundings are also found. When rocketsonde and radiosonde 
observations are found to be coincident with an SSM/T sounding, they are quality tested and 
combined to produce a single sounding extending from sea level to the 0. 5-mb level. This 
composite sounding is then used to calculate the data vector. The measured data vector is 
also determined. These data vectors are finally stored in the SSMT’l'SSMTSTATFILE file with 
the residual error data and the number of coincident data vector pairs. The SSMT^SSMT - 
STATFILE file is used also as an input to MICRO-CHECK and successively updated each day 
with new data. MICRO-CHECK also contains provisions for print -plotting coincident SSM/T 
and radiosonde soundings and listing the statistical data. MICRO-CHECK runs to completion 
when all soundings in the UADB are processed. 


4. 5 D-UPDATE 

Periodically, depending on the number of coincident measurements found in MICRO- 
CHECK, the D-UPDATE program is executed. It runs automatically with input data from 
SSMT*SSMTSTATFILE. Using the data input from this file, D-UPDATE tests the data for 
sample-size representativeness and distribution in temperature. If these tests are satisfac- 
torily passed, an OPTION 1 update may be executed. New D-matrices are generated, tested 
in PROCESSOR-B and used operationally if warranted. After five successive OPTION 1 
updates, tests of OPTION 2 may be performed to determine if the sensor bias and random 
errors exceed the threshold values. When they exceed threshold values, an OPTION 2 update 
is executed by modifying the correlation matrices, testing them in PROCESSOR-B, and 
replacing operational D-matrices with the newly developed ones. 

5. RESULTS 

The SSM/T software has been implemented, tested and verified in the AFGWC UNIVAC 
1110 system using simulated data in the input file PROD*AEROJETMDATA. These data were 
computed from an independent set of soundings selected from the World Data Center A and 
Eastern Meridional high-altitude meteorological soundings. Data were selected for all sta- 
tions in these networks and stratified by latitude and season to complete a nominal DMSP orbit. 
A total of 192 soundings were used in the simulation with 1344 views. Each view was assigned 
latitude and longitude coordinates, surface height, surface temperature and emissivity value 
(over land) or salinity value (in parts per thousand) for views located over the ocean. Antenna 
temperatures were computed for each channel frequency and input to a simulator program 
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with sensor and system data to generate the input file for test and verification purposes. 
Several versions of this file were generated by the simulation program to test and verify 
various system mode characteristics. The inverted soundings produced in the normal proc- 
essing mode were compared to original source reference data by a separate subroutine for 2 
seasons and 6 geographic zones. These included the arctic, mid-latitude, and tropic spring 
in the Northern Hemisphere, and the antarctic, mid-latitude, and tropic fall in the Southern 
Hemisphere. 

The results, obtained from comparing 224 inverted soundings with the source reference 
soundings in each geographic zone, are shown in Figure 4. The standard deviation of the 224 
soundings in each zone-season are also shown in this figure. As shown in Table I, the SSM/T 
sensor can measure brightness temperatures averaged over the weighting functions shown in 
Figure 2 with an rms accuracy of 0. 4 to 0. 6 K for integration time of 2. 7 seconds. The 
results shown in Figure 4 indicate the SSM/T rms uncertainty averaged over the 15 mandatory 
levels is of the order of 2 K. Larger uncertainties occur at the 1000-mb heights since the 
inferred soundings are biased to the 1000-mb forecast height. A qualitative comparison of 
the rms uncertainty in the inferred soundings with the standard deviation of the source sound- 
ings reveals that soundings will be obtained from the SSM/T measurements with an accuracy 
considerably better than estinnates based on climatology. 


TABLE I. CHANNEL PARAMETER DESIGN SPECIFICATIONS 




F requency 

Bandwidth 

NETD 

Channel 

Polarization 

(GHz) 

(MHz) 

(OK) 

1 

) 

50. 5 

400 

0. 6 

2 


Principally 
' Horizontal 

53. 2 

54. 35 

400 

400 

0. 4 
0. 4 

3 


4 

) 

54. 9 

400 

0. 4 

5 



58. 4 

115 

0. 5 

6 


Orthogonal 
to Channels 

58, 825 

400 

0. 4 

7 

J 

1-4 

59. 4 

250 

0. 4 


TABLE II. SCAN PARAMETERS 


Scan Type 

Cross-Track Positions 

Calibration Positions 

Instantaneous Field of View 

Total Cross-Track Scan 

Total Scan Period 

Integration Time (Cross -Track and Calibration Positions) , . 
Sync Mode 


Cross-Track Nadir 
7 

2-Cold Space and 300°K 
14° 

36° 

32 Seconds 

2. 7 Seconds 

Auto Mode or On-Sync 
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FIGURE !• Scan Geometry 
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FIGURE 2. Weighting Function, 
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FIGURE 3. SSM/T Software Diagram 
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